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ABSTRACT 



We report the possible detection (99.3% of statistical significance) of 
redsliifted Fe iron line emission in the X-ray afterglow of Gamma-ray burst 
GRB970508 observed by BeppoSAX. Its energy is consistent with the redshift 
of the putative host galaxy determined from optical spectroscopy. The line 
disappeared ~ 1 day after the burst. We have also analyzed the spectral 
variability during the outburst event that characterizes the X-ray afterglow of 
this GRB. The spectrum gets harder during the fiare, turning to steep when 
the fiux decreases. The variability, intensity and width of the line indicate that 
the emitting region should have a mass ^ O.5M0 (assuming the iron abundance 
similar to its solar value), a size of ~ 3 x lO^^cm, is distributed anisotropically, 
and is moving with sub-relativistic speed. In contrast to the fairly clean 
environment expected in the merging of two neutron stars, the observed line 
properties would imply that the site of the burst is embedded in a large mass 
of material, consistent with pre-explosion ejecta of a very massive star. This 
material could be related with the outburst observed in the afterglow 1 day after 
the GRB and with the spectral variations measured during this phase. 

Subject headings: gamma rays: bursts - X-rays: general - line: formation 
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1. Introduction 

Distance - scale determination of Gamma-ray bursts (GRB) has been one of the most 
important achievements of astrophysics in recent years. Accurate and fast locahzation of 
the prompt and afterglow emission bursts (i.e. |Costa et al. 1997| ) by the Italian/Dutch X-ray 



satellite BeppoSAX (|Piro, Scarsi fc Butler 1995| , Poella et al. 199^ ) led to the identification 



of optical counterparts ( [van Paradijs et al. 1997|) and ultimately to spectral measurements 



of a redshift ( [IVletzger et al. 1997| ). While the extragalactic origin of GRB has gathered 



solid evidence in its support, the source of the large energy implied by their distance is still 
speculative. 

Direct information on the nature of the central engine of GRB can be derived by 
studying the nearby environment; for example using line spectroscopy. The star forming 
region and massive pre-explosion winds associated with the hypernova scenario (Paczynsky 



1998|) imply a mass rich environment down to short (~ IQ^^cm) distance scales. In contrast, 



NS-NS merging should happen in a fairly clean environment, because such objects are 
expected to form with significant speeds which can lead them at least a few parsecs away 
from the regions of star formation in which they formed. Current line measurements are 
inconclusive, because all the spectral features observed so far are signatures produced 
by the host galaxy rather than at the burst site. GRB980425 is an exception, but its 
tantalizing association with SN 1998bw ([Galama et al. 1998|, [Kulkarni et al. 1998|) requires 



clarification. 

The measurement of X-ray lines emitted directly by the GRB or its afterglow could 
provide a direct measurement of the distance and probe into the nature of the central 
environment ( |Perna fc Loeb 1998| ). Fe K line emission is very promising in this respect. 



Theoretical computations have been carried out by several authors in the framework of 
GRB ( llVleszaros fc Rees 1998| , [Boettcher et al. 1999| , |Ghisellini et al. 1998| ). Most of this 
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work deals with features generated within the interstellar medium (ISM) of the host galaxy, 
i.e. on scale larger than several parsecs. However, even in the favourable case of dense 
regions of stellar formation (e.g. the site of GRB formation in the hypernova scenario), 
iron emission lines would be hard to detect with current X-ray instrumentation. Lately, 
Meszaros and Rees (1998) have shown that the circumburst environment created by the 
stellar wind before the explosion of the hypernova could yield a line of substantial intensity. 
A similarly favourable situation should be expected in a variation of the hypernova 
scenario, - the SupraNova ( |Vietri fc Stella 1998| ), where the GRB is shortly preceded by a 



supernova explosion with the ejection of an iron-rich massive shell. In these cases, and for 
a isotropic distribuition of material around the source, a tight upper limit on the column 
density N^ < lO^'^cm"^ derives from the requirement that the Thomson optical depth be 
less than a few, to avoid smearing out the the short time-scale structure of the GRB by 
Thomson scattering ( Poettcher et al. 1999|) . Furthermore, Fe absorption features of stregth 
comparable to that in emission should also be detected. 

It is also conceivable that the impact of the relativistic shell that produced the original 
GRB on these ejecta could provide an additional energy input in the afterglow. 

GRB970508 is the most promising candidate of the BeppoSAX afterglows for the 
search of line emission. It is characterized by a large outbursting event during its afterglow 
phase (|Piro et al. 1998| ) and has the highest signal to noise ratio of the BeppoSAX GRB 



afterglows. In this Letter we present results of a spectral analysis of the afterglow of 
GRB970508 and discuss their implications on the nature of the central engine of GRB. 
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2. Observations of the X-ray afterglow spectrum of GB970508 

The Gamma-ray burst GRB970508 was detected by the Gamma-ray Burst Monitor 
(GRBM) aboard BeppoSAX ( piro et al. 1998| ). The Wide Field Camera (WFC) provided 



a few arcmin position acquired with the Narrow Field Instruments 6 hours after the burst, 
which led to the identification of the X-ray afterglow associated with the GRB. Three 
further observations were carried out, the last performed 6 days after the burst (Fig.l, 
Tab.l). 

LECS(0.1-3.5 keV) and MECS(1. 8-10.0 keV) spectra of the observations were obtained 
using standard procedures (BeppoSAX web pages: www.sdc.asi.it), with 15-20 counts 
(source -|- background) per bin. 

Inspection of the light curve (Fig.l) shows that in the first half period of observation 
number 1 the fiux decays, followed by a rising trend in the second half. We have thus split 
this observation into two parts (la & lb) at the time of the minimum of the light curve, to 
investigate spectral evolution in different states, producing 5 spectra in total. 

2.1. The iron K^ line 

We fitted each spectrum with a power law with variable absorption, keeping the relative 
normalization of LEGS vs MEGS free to vary in the range 0.6-0.9. This model provides a 
good fit to all but dataset la. In this case we have xl = 2.4 for u = 8 degrees of freedom, 
which corresponds to a probability of 1.4% [| 



^Note that GB970508 has the highest S/N ratio of the BeppoSAX GRB afterglows, so 
that the number of spectra with comparable statistical weight is limited to those of data set 
1. We also remind that the afterglow of GB970228 had an higher fiux, but the S/N ratio 
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However the major contribution to x^ is not randomly distributed but is concentrated 
in a line feature around 3.5 keV (Fig. 2a; note that in figure we show for clarity only data 
points below 2 KeV for LEGS and above 2 keV for MEGS). 

The possible association of this feature with Fe K line emission is obvious, considering 
that its energy coincides with Fe I-Fe XX Ka emission at the redshift z=0.835, i.e. that 
of the putative host galaxy (Bloom et al. 1998). We have therefore added to the power 
law model a narrow line (width < 1.5 x FWHM~ 0.5 keV) at the energy of 3.5 keV. 
The relative normalization LEGS/MEGS - that depends primarily on the position of the 
source in the LEGS detector has, in all datasets, an error greater than the nominal range 
0.6-0.9. We have hereafter fixed it to 0.8, the value averaged over all data sets. The best 
fit now yields xl = 0-93 with u = 7. The improvement of Ax^/xl = 13.6 corresponds to a 
confidence level of 99.3% (F test). We stress that this confidence level relies on the a priori 
knowledge of the energy of the Fe line available from an independent redshift measurement . 
Were it not available, the confidence level for the presence of a line with free energy would 
have been 97%. 

The best fit line intensity is Ipe = (5 ± 2) x 10~^ ph cm~^s^^ (hereafter all errors and 
upper limits correspond to 90% confidence level for a single parameter). Upper limits on 
line intensity in datasets lb (Fig. 2b), 2 3 and 4 are 2.5, 2, 2.5 and 0.5 xlO~^ ph cm~^s~^ 
respectively. 

The power law parameters obtained in dataset la and lb are consistent with each other 
(see Tab. 1) and the absorption column densities are consistent with the hydrogen column 
in the line of sight of our Galaxy, i.e. the intrinsic absorption is consistent with zero. We 
next performed a simultaneous fit to both datasets with the same continuum parameters. 



was lower, because data of one MEGS only were usable( posta et al. 1997|) 
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and the addition of a line to dataset la only. In this case we left the energy of the line 
as free parameter, deriving E = 3.4 ± 0.3, corresponding to (6.2 ± 0.6) in the source rest 
frame. The error range is too large to identify the ionization state of iron. In fig. 2 (inset) 
we show the contour plot of I^e vs E. 

We note that the absence of the line in dataset lb explains why a power law model 
provides an acceptable fit to the summed spectrum ( Pwens et al. 1998|) . We have also 
investigated whether the spectral continuum could be fitted by a black body or a thermal 
bremsstrahlung. We find that thermal models provide worse results than that of a power 
law continuum, as already shown in the case of GRB970228 ( p^'rontera et al. 199(j| ). 

2.2. Spectral variability in the afterglow outburst 

Spectral variability of the afterglow spectra supports the idea that a refreshed energy 
input produced the outburst observed in the light curve. The spectral index in the first 
observation (Tab. 1) is a = 1.5 ± 0.6. In the second and third observation the flux lies 
above an extrapolation of the light curve expected from a power law evolution from the 
primary event. The spectrum is harder (a = 0.4 ± 0.6), indicating electron reacceleration, 
as expected in a newly formed shock. In the last observation the flux decreases and the 
spectrum is very steep (a = 2.2 ± 0.7): the new energy source has been exhausted. We 
have noticed previously that spectra la and lb are consistent with each other. Since 
the transition la-lb identifies the beginning of the outburst, one would then expect 
an hardening of the spectrum in lb. However such an effect could be masked by the 
continuation of the power law decay law (Fig.l), which accounts for about 60% of the flux 
observed in dataset lb. 



3. Origin of the line and constraints on the emitting region 

The minimuiii amount of mass Mmm present can be estimated (Lazzati et al. 1998) 
simply by counting how many ionizing photons a given iron atom can absorb, assuming 
that recombination is much faster than the time interval between successive ionizations, 
and by demanding that the existing iron atoms thus account for the fiuence in the line: 

ylO ^ ph s ^ cm ^ j \\P s V10^° cm) q 
where Ap^ is the iron abundance in units of the solar value, T is the line duration, 
R = 10^^-Ri6 cm is the distance of the medium from the central source, and q < 0.1 is the 
fraction of the total energy E = IO^^-E'52 erg which is recycled into the line. Since the line 
disappears after about 10^ s, we expect the emitting region to have a size D ^ 3 x 10^^ cm; 
in fact, since the line width is mostly instrumental, we can bar the existence of relativistic 
effects. Then the typical density is n = 5 x 10^ AplRfg cm.~^, and typical Thompson and 
iron optical depths are tt = 10 AplRfg and Tp^ = S-Rfg. 

/^From the above we deduce that the emitting material cannot lie along the line of 
sight: for Riq > 1, tt ^ 1, and Thomson scattering would smear out the time structure of 
the burst, contrary to observations ( [Boettcher et al. 1999| ). Even if Ape ^ 1, we would still 
have Tpe ^ 1, which implies an iron edge not seen in the spectrum. For Riq < 1, we notice 
that the fireball would load itself with so many baryons to spoil its relativistic expansion, 
again contradicting observations. 

Having established that the emitting medium lies sideways from the line of sight, we 
also notice that we cannot have R ^ D, because then the line would appear much later 
than 1 day. We conclude that R ^ D. 

The ionization parameter, using the minimum density derived above, is found to be 

e = -^ = lO^^Fe .^so ^ -1 ^1^' ^ 10' Ap,t-'-'R^,' (2) 

nR'^ 10^^ erg s ^ 
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where L is the luminosity in the range 0.016-10 keV, whose temporal evolution is 
approximately described by a power law ( piro et al. 1998]) . /^From Eq. (2) one sees that 



for early times (t ^ 10^ s), the material will be fully ionized, and thus the line may be 
due to recombination, while at later times, when ^ ^ 10^, fluorescence may account for 
the emission ([Kallmann fc McCray 19^ , [Hirano et al. 1987|) . The exact moment of the 



transition cannot be determined, because we established above that the emitting material 
lies sideways with respect to the observer, and we cannot exclude some amount of beaming, 
but the following discussion still stands. At early times, the emitting material will probably 
reach close to the Inverse-Compton temperature in the burst radiation fleld, which clearly 
exceeds 10^ °K; for these high temperatures and low densities, recombination is not an 
effective emitter of line photons, and it is most likely that most of the line we see is emitted 
by fluorescence in the time range 10^ s < t < 10^ s, during which the total observed fluence 
is ~ 10^^ erg. Putting this into Eq. |I|, we flnd a minimum mass of 

M^,„ = 0.5MqA]^i (^) , (3) 

which is our best estimate for the total mass present. 

3.1. The site of the line emission 

In principle, the line might arise from the normal interstellar medium surrounding the 
burst volume, or from a normal stellar wind emitted by the burst progenitor. However, 
a very large density {n ^ 10^ cm"^) is necessary to account for the observed spectral 
feature; this is obvious when one considers that the line reported here is as bright as a 
typical one from clusters of galaxies at z ~ 0.1 <C 2^970508 • Winds of massive stars are 
also unlikely to be responsible for the iron line: even for rather extreme wind parameters 
(M = m_4 X lO^^M© yr"^, and f = f 2 x 100 km s^^), the wind density exceeds the minimum 
density only for Riq < 0.1(m_4A^e/'^2)^^^- However, for this small distance, the ionization 
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parameter (Eq. ^ is ;^ 10^ even a day after the burst and, even including beaming, the 
wind material will be totally ionized; under these conditions, only recombination can 
provide the line photons, but the recombination rates are very small {tree > 10^ s), and thus 
the conditions of instantaneous de-excitation under which Eq. |l] was established are no 
longer satisfied: much larger masses are required than can be accounted for by the winds. 
Furthermore, such small distances {Riq < 0.1) cannot explain the line duration. We thus 
conclude that the material surrounding the burst must have been pre-ejected. 

Ape ?^ 1 is likely to be close to an upper limit to the iron abundance of the ISM 
or of a stellar wind at 2; ~ 1. Scenarios concerning mergers of binaries then begin to 
appear unlikely, regardless of whether the binary is made of two neutron stars, a black 
hole and a neutron star, or a black hole and a white dwarf, because it has been established 
(Meszaros and Rees 1992) that they cannot pre-eject, even under the rosiest assumptions, 
more than si 10~^ Mq due to their tidal interactions, about three orders of magnitude less 
than the limit in Eq. |I]. We remark that the SupraNova (|Vietri fc Stella 1998|) provides 



the most appealing scenario, because of the large amount of iron material ejected by the 
supernova explosion preceding the burst. In fact, theoretical computations of the explosion 
of massive stars indicate that ~ O.IM© of iron are ejected by the SN ( Woosley fc Weaver 



1995|, [I'hieleman, Momoto fc Hashimoto 1996|). In the case of SN1987a, whose progenitor 



had M Ri 2OM0, a similar amount of iron was derived from optical observations ( panzige"r 
fc Bouchet 19931) . 



In future, it is conceivable that simultaneous observations of the Ly-a line may help to 
determine Ape, and thus uniquely determine the total illuminated mass, and by implication 
the adequacy of different formation scenarios. 
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4. Summary 

In this paper we have obtained the following results: 

• We have searched the X-ray spectrum of GB970508's afterglow for an iron line, 
located at the system's redhift (z = 0.835); we found such a line with limited statistical 
significance (99.3%) in the early part (first 16 h) of the afterglow; the line decreases 
in the later part of the observations (~ 1 day after the burst) by at least a factor 2, 
enough to make it undetectable with current apparatus; 

• Simultaneously with the line disappearance, the X-ray flux both rises and hardens 
{a = 0.4 ±0.6, while a = 1.5 ±0.6 before the reburst), consistent with the appearance 
of a new shock. Then, at the end of the outburst, the spectrum steepens. 

We showed that these observations indicate that a mass ~ 0.5MqApI should be 
located at a distance of 3 x 10^^ cm and sideways with respect to the observer, and is 
moving at subrelativistic speed. In order to reach such distance, this material must have 
been pre-ejected by the source originating the burst, shortly (perhaps a year, for typical 
SN expansion speeds) before the burst. We stress that these observations contain two 
coincidences: on the one hand, this is the only burst in which a reburst and a line have 
been observed by BeppoSAX; on the other, the iron line disappears exactly at the moment 
of the reburst. We also point out that a line feature, with a similar significance, has been 
found by ASCA in another burst, GRB970828, which also shows an event of rebursting 
during the X-ray afterglow ([Yoshida et al. 19981) . 



We hope the report of this result will help shape future observations by X-ray satellites 
with superior observing capabilities, such as XMM, AXAF, ASTRO-E. 
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Fig. 1.— The X-ray (2-10 keV) light curve of the X-ray afterglow of GRB980508. The 
dashed line represents the extrapolation from the X-ray data points of the burst (Piro et 
al. 1998). For the spectral analysis data from the first pointing was divided in two parts, 
corresponding to the transition between the decay and turn up from the light curve, creating 
in total 5 spectra (la,lb-4). 

Fig. 2. — The spectra (in detector counts) of the afterglow of GRB970508 from data set la 
(panel a) and data set lb (panel b). The continuous line represents the best fit power law 
continuum in the range 0.1-10 keV, excluding the interval 3-4 keV. For clarity we show in 
figure LEGS data below 2 keV and MEGS data above 2 keV. The inset in panel a) shows 
the contour plot of of the line intensity vs energy. Gontours correspond to 68%, 90% & 99% 
confidence levels for two interesting parameters. 
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Table 1. Log of Observations and fit with a power law 
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